Red Nuggets at z ~ 1.5: Compact passive galaxies and the 
formation of the Kormendy Relation 
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ABSTRACT 

We present the results of NICMOS imaging of a sample of 19 high mass passively evolving 
galaxies -with 1.2 < z < 2, taken primarily from the Gemini Deep Deep Survey (GDDS). Around 
80% of galaxies in our GDDS sample have spectra dominated by stars -with ages > 1 Gyr. Our 
rest-frame i?-band images sho-w that most of these objects have compact regular morphologies 
-which foUo-w the classical R^/^ la-w. These galaxies scatter along a tight sequence in the size 
vs. surface brightness parameter space -which defines the Kormendy relation. Around one-third 
(3/10) of the massive red objects in the GDDS sample are extraordinarily compact, -with effective 
radii under one kiloparsec. Our NICMOS observations allo-w the detection of such systems more 
robustly than is possible -with optical (rest-frame UV) data, and -while similar systems have been 
seen at 2: > 2, this is the first time such systems have been detected in a rest- frame optical 
survey at 1.2 < z < 2. We refer to these compact galaxies as 'red nuggets', and note that 
similarly compact massive galaxies are completely absent in the nearby Universe. We introduce 
a ne-w 'stellar mass Kormendy relation' (stellar mass density vs size) -which we use to single out 
the effects of size evolution from those of luminosity and color evolution in stellar populations. 
The 1 < z < 2 passive galaxies have mass densities that are an order of magnitude larger then 
early type galaxies today and are comparable to the compact distant red galaxies at 2 < z < 3. 
We briefly consider mechanisms for size evolution in contemporary models focusing on equal- 
mass mergers and adiabatic expansion driven by stellar mass loss. Neither of these mechanisms 
appears able to transform the high-rcdshift Kormendy relation into its local counterpart, leaving 
the origin and fate of these compact 'red nuggets' unresolved. 

Subject headings: galaxies: elliptical, galaxies: fundamental parameters, galaxies: evolution 



t Department of Astronomy & Astrophysics, University 
of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4 

■•■Observatories of the Carnegie Institution of Washing- 
ton, 813 Santa Barbara Street, Pasadena, CA 91101 

^Centre for Astrophysics and Supercomputing, Swin- 
burne University of Technology, 1 Alfred St, Hawthorn, 
Victoria 3122, Australia 

^Max-Planck-Institut fiir extraterrestrische Physik, 
Garching, Germany 

°Herzberg Institute of Astrophysics, National Research 
Council, 5071 West Saanich Road, Victoria, British 
Columbia, V9E 2E7, Canada. 

® Department of Astronomy /Steward Observatory, Uni- 



versity of Arizona, 933 N Cherry Ave., Rm. N204, Tucson 
AZ 85721-0065 

*DSM/ DAPNIA/ Service d'Astrophysique, CEA/ 
SACLAY, 91191 Gif-sur-Yvette Cedex, France 

^Gemini Observatory, Hilo, HI 96720 

^The Department of Astonomy and Astrophysics, Uni- 
versity of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637 

*Dept. of Physics and Astronomy, San Francisco State 
University, 1600 HoUoway Avenue, San Francisco, CA 
94132 



1. Introduction 

The formation mechanism of eUiptical galaxies 
has long been controversial and remains a key test 
of more general galaxy formation models. The 
origin al 'nature' (JEggen. Lvnden-Bell fc Sandaed 
([1962) monolithic collapse) vs. 'nurture ' (for- 



mation through mergers, fe.g.. ISchweizer 111987 : 
Searle fc Zinnl(l978t [Tbomre k Toomrdll972l ) de- 



bate is still with us, but is now set in a ACDM 
cosmological context which attempts to connect 
the stellar component of galaxies to an underly- 
ing evolutionary picture for the clustering of dark 
matter halos. Testing this model requires study- 
ing the evolution of galaxies over a large redshift 
range. 

A wide range of selection techniques have been 
effective in selecting galaxies in various redshift 
ranges on the basis of their current star for- 
mation rates (e.g. Lyman break galaxies, sub- 
mm sources etc), or from the spectral signa- 
tures of passively evolving old stellar populations 
(e.g., extremely red objects (EROs) and other 
color selections). The most massive local ellip- 
tical galaxies have the oldest s tellar populations 
([Gallagher. Hunter fc Tutukovl Il984l ). so identi- 
fying the progenitors of local early-type galaxies 
within the high-redshift galaxy population is of 
particular interest. There is a consensus that 
the mass density in the red sequence is evolv- 
ing strongly in the 1 < z < 2 r ange (GDDS 
Paper VIII 
per III 



Abraham et "ail 120071: GDDS Pa- 



Glazebrook et all 2004; 



Fontana et al.l 



2004t iRudnick et al.ll2003l ). a process that contin 



ues a t redshifts below unity as well (jFaber et al. 



20071 : iBellet al.M2004f) . altho ugh the magnitude 
of the evolution i s uncertain ( Brown et al.ll2007t 
20031 1 . Massive morphologically- 



Chen et al. 



confirmed elliptical galaxie s have been found up 
to z = 2 (GDDS Pa per IV, iMcGarthv et~alll2004l : 
Cimatti et al.l 120041 ) with spectra consistent with 
formation epochs up to z > 5. These observations 
were in in direct contradiction with early ACDM 
models where stellar mass assembly traced the 
build up of cold dark matter haloes, although 
additional feedback mechanisms on the baryons 
have more r ecently been able to bet ter account 
for this (e.g.. lKang. Jing. fc Silkll2006h . A compli- 
cation recently added to this picture is the obser- 
vation that the space density of ellipticals is found 



to evolve strongly over 1 < z < 2 (jPaper Vllj ) 
even while their stellar populations evolve weakly, 
suggesting that one must be careful to decou- 
ple morphological evolution from evolution of the 
underlying stellar populations. This is seen at 
higher redshifts also, where the paucity of pas- 
sively evolving galaxies at z > 2 in deep J — K 
and 3.5 nva selected samples ( Kriek et al.l l2006t 
Labbe et alJ 2005; C imatti et al.ll2002f ) shows that 
the assembly epoch for the red sequence may be 
decoupled from the epoch of the earliest star for- 
mation. Studies of star formation history and 
morphology can only go so far in unraveling the 
puzzle of galaxy formation; dynamical and chemi- 
cal probes are needed to connect progenitors to de- 
scendants. Clustering signatures offer one dynam- 
ical approach to connecting progenitors to descen- 
dants an d the strong clustering of t h e passive red 
galaxies (iDaddi et all l2005aL 120041 : iBrown et al 



l2003t iMcCarthv et al.ll200ll ) stronglv suggest that 
they are linked to today's massive ellipitical galax- 
ies. 

Theoretical attempts to explain these observa- 
tions have resulted in greatly improved ACDM 
models which decouple mass assembly from this 
stellar population downsizi ng. An example is 
the s emi-analytic model of iDe Lucia fc Blaizot 
(J2007[ ). Here the small ellipticals and their 
stars form early by disc mergers. Massive el- 
lipticals can then grow bigger and more numer- 
ous at late times through dissipationless or dry 
merging. This may even have been observed 
(JBell et al.l l2006a) though there is some disagree- 
ment as to whether the ACDM merger rate is 
high enough (JBundv et al.ll2007[ ). At this stage 
it is perhaps fair to say that dry merging sim- 
ulations show that i t does not disrupt ellipti- 
cal scaling relati ons (IBov lan-Kolchin et al. 120061 
2005; Gonzalcz- Garcia fc v an Albada 2003) as one 
might naively expect. However only a limited 
number of simulations of this process have been 
done and they have not yet been incorporated 
into cosmological models in a detailed way such 
that they can be compared with data (e.g., num- 
bers, sizes and masses of galaxies). Further it is 
not clear that a dry merging hierarchy consistent 
with cosmological downsizing can also be made 
consi stent with the evolving mass -metallicity rela- 
tion (jPipino and Matteuccil 2008 ). A co ntrasting 
picture is painted by iNaab et al.l (J2007[ ) using a 



SPH model of individual systems. They argue 
for a formation mode dominated by something 
very close to early monolithic collapse, but in a 
ACDM cosmological context, with mergers (along 
with accretion) playing only a minor role in stellar 
mass growth at late times. 

High spatial resolution studies of the morpholo- 
gies and structures of passive galaxies offer one 
approach to gauging the importance of recent ma- 
jor merger events. A number of studies with the 
Hubble Space Telescope (HST) have shown that 
half or more of red galaxies in color-selected sam- 
ples have simple early type morphologies. Most of 
these studies are confined to redshifts of ~ 1.5 and 
less, an d the early- type fraction varies from ~ 50% 
to 70 % (JMoustakas et al.ll2004l: lYan fc Thompson 
2003f ). At higher redshifts a significant frac- 
tion of the red galaxies appear to be discs (e.g.. 
Paper VII]| . iFontana et al.l 12004 ). Understanding 
the connection between these two classes of ob- 
jects naturally focuses on the importance of merg- 
ers, since nearly equal-mass mergers are thought 
to transform discs into spheroids. Mergers, both 
gas-rich and dissipationless, are also thought to be 
important in the growth of the red sequence and 
evidence, both direct and indirect, supports that 
this i s occuring at int ermediate and low redshifts 
(e.g.. lBell et al.l2006bL and the references therein). 
It appears that much of the high-redshift merg- 
ing activity may be of the dissipationless variety 
where the main effect of merging is to reorga- 
nize existing stellar population without triggering 
new star formation. It is difficult to envision how 
this might operate unless the merging systems 
are thems elves gas-poo r, which is not generally 
expected ( van DokkumI 12005). In any case, the 
signatures of such 'dry' mergers are difficult to 
detect at high redshifts. 

Recently, several imaging studies have shown 
that red galaxies at z > 1 appear smaller than 
their likely p resent-day descendants with the same 



stella r mass (jLonghetti et al.l l2007t iCimatti et al. 
20081 ). The implications of these observations 



are seen most clearly in the structural and dy- 
namical scaling relations, the Fundamental Plane 
and its projections (the Faber-Jackson (1976) and 
Kormendy (1977) relations). In the present pa- 
per we explore the nature of the Kormendy re- 
lation, (mean surface brightness within the effec- 
tive radius, (/i)e, versus effective radius. Re). This 



is the most observationally accessible projection 
of the fundamental plane at high-redshift. Our 
analysis spans the redshift range 1.2 < z < 2 
using HST NICMOS observations of a sample 
of quiescent high-redshift galaxies taken mainly 
from the Gemini Deep Deep Survey (GDDS Pa- 
per I, |Abraliam_eialJ|200j). We present NICMOS 
F160W images for ten of the twe nty z > 1.3 pas- 
sive red galaxies from IPaper IVl . These systems 
all have spectra dominated by old stellar popula- 
tions. This extends to higher re dshifts (z > 1.7) 



than t he earlier NICMOS work of iLonghetti et al. 
(|2007l ) from the Munich Near-IR Cluster Sur- 



vey (MUNICS, iDrorvet al.ll200lf l. We also in- 
de pendently analyze the archival NICMOS data 
of ILonghetti et al.l ( 20071 ) in the redshift range 
1.2 < z < 1.7 to supplement our sample and con- 
firm their findings. At the higher redshifts pre- 
vious findings of compact galaxies were based on 
optical data obtained with the Adv anced Camera 



for S urveys (ACS) onboard HST (jCimatti et al. 
I2OO8I ). Our use of NICMOS allows us to more ro- 
bustly show that the old components in the galax- 
ies are truly compact. Finally, we are able to unify 
the optical and infrared work by introducing a new 
'stellar mass Kormendy relation' which we use to 
better quantify evolution in the sizes of early-type 
galaxies as a function of stellar mass over the red- 
shift range 1 < 2 < 2. We briefly examine the 
likelihood that dry mergers explain such size evo- 
lution, and examine whether an alternative pro- 
cess, adiabatic expansion, might be more impor- 
tant. We describe the observations in section [2l 
our analysis in section [Sj and present our results 
in section H) In section [5] we discuss the impli- 
cations of our observations for simple models for 
galaxy size growth. Throughout we use standard 
cosmological parameters; Hq = 70 km s~^ Mpc~^, 
f^m = 0.3, r^A — 0.7. Unless stated otherwise, all 
magnitudes are based on the AB system. 

2. Description of the Observations 
2.1. Sample definition 

Our sample of galaxies was taken mainly from 
the GDDS, crafted to sample the galaxy popula- 
tion in the critical 1 < z < 2 in terval with an 
emphasis on red galaxies ( Paper j ) . While mod- 
est in area (120 square arcminutes), the survey is 
spread over four independent and representative 



sightlines. Redshifts for ~ 300 galaxies brighter 
than J(Vega) = 24.5 were obtained from 30-hour 
long integrations using the GMOS spectrometer 
on Gemini North. This magnitude limit corre- 
sponds to the stellar mass of 2.5 x 10^° Mq for 
a galaxy with the redshift of formation Zf — 10 
and maximally ol d stellar po pulation observed at 
redshift z = 1.5 (jPaper IIll ) . We classified the 
galaxies on the basis of their spectra, depending 
on whether they were dominated by active star for- 
mation, stars older than ^ 1 Gyr, intermediate age 
(0.3 — 1 Gyr) populations, or a mix of these types. 
Of the 302 galaxies with redshifts, 47 have spec- 
tra dominated by old stars, and twenty of these lie 
at redshifts beyond 1.3. Spectra of these twenty 
galaxies and estimates of their ages and formation 
redshifts are presented in IPaper IVl . Deep J-band 
images of the GDDS galaxies at z < 1.7 with the 
ACS on HST reveal that the correlation between 
spectral type, and hence stellar content, and mor- 
phological class seen at present is strong at these 
redshifts. Nearly all of the GDDS galaxies with 
passive spectral classes have compact morpholo- 
gies consistent with early Hubble types, while the 
actively star forming galaxies have a morphologies 
that range from simple disks to complex struc- 
tures indicative of ongoing mergers. The GDDS 
galaxies discussed in this paper are a subset of the 
GDDS galaxies having spectra dominated by old 
stars (class "001" from lPaper J ) and z > 1.3. The 
key properties of this sample are given in Table [1] 

Our primary sample of ten galaxies is drawn 
from the GDDS and determined by the number 
of available orbits and the desired depth of NIC- 
MOS imaging. The targets were selected ran- 
domly, with the exception of the two (12-5869 
and 12-5592) that could be covered in a single 
pointing. We also analyzed archival data from 
the MUNICS survey for nine additional galaxies 
with properties similar to th ose of our GDDS sam- 
ple. iLonghetti et al.l (20051) analyzed spectropho- 
tometric data set for these galaxies from the ncar- 
infared spectroscopic follow-up of a complete sam- 
ple of bright {K < 18.5) EROs (R - K > 5.3) 
selected from the MUNICS survejtJ. Low resolu- 
tion spectroscopic and photometric data revealed 
stellar masses greater than 10^^ M0 and dominant 



old stellar population for all objects in the sample 
(see Table [5]). As will be described below, this 
additional data provided us with a useful check of 
our methodology by allowing us to compare results 
fr om our analysis pipelin e against those published 



Longhetti et all (|2007l ) 



2.2. NICMOS Observations 

The ten GDDS galaxies were observed with 
Camera 3 on NICMOS using the F160W filter. 
Each individual exposure was 896 seconds in du- 
ration with multiple samples using the STEP64 
read pattern. A single orbit contained three expo- 
sures and we observed each target over four HST 
orbits for a total integration time of 10740 sec- 
onds. Two of the fields overlapped and the im- 
ages for targets 12-5869 and 12-5592 have twice 
the exposure time of the others. These objec ts are 
discussed in detail in iMcCarthv et al.l (|2007n . We 
dithered in non-integer pixel steps between each 
exposure. The individual frames were dark cor- 
rected, sky subtract ed and combined usiii g the 
DRIZZLE algorithm (JFruchter fc Hookl2002h with 
a final pixel size of 0'.'12. Residual sky levels in the 
final mosaics were derived from Gaussian fits to a 
histogram of sky values and were subtracted. 

As noted above, we also re-analyised nine galax- 
ies from the MUNICS sam ple of red galaxies de- 
scribed in ILonghetti et all i|2007l ). The MUNICS 
data set was obtained using Camera 2 on NIC- 
MOS, and is thus more finely sampled, and some- 
what shallower, than our NIC3 images. As de- 
scribed below, analyzing this NIC2 data allowed us 
to explore, and ultimately rule out, the possibility 
that the coarser sampling of our NIC3 data might 
lead to poor model fits and spurious sizes. We 
retrieved the pipeline-processed individual NIC2 
images from the HST archive. We then corrected 
each image for residual pedestal effects and com- 
bined them into mosaics using the DRIZZLE al- 
gorithm with a final pixel size of 0'.'05. The prop- 
erties of the nine galaxies in this sample are sum- 
marized in Table 2. 

3. Analysis 

3.1. Surface brightness profiles 



^This is actually a blank field survey, the intention was to 
find high-2 clusters from deep wide-field near-IR imaging. 



Using the Galf it software package (jPeng et al 



20021 ) ■ we derived two-dimensional (2D) surface 



Table 1 
Properties of the ten galaxies in GDDS sample 



ID 


z 


Mass'' 


Age^ 






[10" Mq] 


[Gyr] 


12-5592 


1.623 


1.16 ± 0.27 


l-lt°ni 


12-5869 


1.51 


3.14 ± 0.43 


1 9+0.6 
^•^-0.2 


12-6072 


1.576 


0.59 ± 0.27 


l.Gtn 


12-8025 


1.397 


1.25 ± 0.39 


o.8t°;? 


12-8895 


1.646 


3.18 ± 0.44 


9 c+0.3 


15-4367 


1.725 


0.56 ±0.15 


2.lt"n1 


15-5005 


1.845 


0.67 ± 0.24 


o.5+t!-I 


15-7543 


1.801 


1.06 ± 0.30 


0.9+°;^2 


22-0189 


1.49 


2.85 ± 0.98 


o n+O.T 
O.U_o.2 


22-1983 


1.488 


1.34 ± 0.53 


1 1+3-^ 


•^ GDDS 


mass 


estimates are 


based on 


the Baldrv & Glazebrookl (|2003|) 


IMF, and 


taken from 


Paper 


la 













Minimum galaxy ages from lPaper Y\ 



Table 2 
Properties of six massive galaxies in MUNICS sample'' 



ID 



Mass'' 
[10" M© 



Age 

[Gyr] 



S2F5.109 


1.22 


5.94 ± 0.95 


1.7 ±0.3 


S7F5.254 


1.22 


4.68 ± 0.16 


5.0 ±0.1 


S2F1.357 


1.34 


4.65 ± 0.40 


4.0 ± 0.1 


S2F1.389 


1.40 


2.15 ±0.86 


3.0 ± 0.5 


S2F1.511 


1.40 


2.07 ± 0.89 


1.3 ±0.3 


S2F1.142 


1.43 


4.06 ± 0.94 


2.2 ± 0.2 


S7F5.045 


1.45 


3.58 ± 1.10 


1.7 ±0.3 


S2F 1.633 


1.45 


3.52 ± 0.51 


4.0 ± 0.5 


S2F 1.443 


1.70 


3.58 ± 1.48 


3.5 ± 0.3 



"" from lLonghetti et all (|2007l ) 

'° MUNICS m a ss est imates are taken from 
Longhetti et "!!! (l2005l Salpeter IMF ) and 



transformed to iBaldrv fc Glazebrookl (|2003|) 
IMF following the relation given in lPaper III 
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Fig. 1. — NIC3 images and the results of our 2D fitting with Galf it for our sample of 10 GDDS galaxies with 
1.3 < z < 2 and spectra dominated by old stars. The three columns present the drizzled F160W image, the 
best fitting R^/^ model, and the residuals. The residual images have been scaled by a factor of 10 compared 
to the data and models to bring out faint features. The bars at the bottom are one arcsecond in length. 
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Fig. 2 . — NIC2 images and the results of our 2D fitting with Galf it of the six galaxies from lLonghetti et al.l 



( 20071 ). The three columns present the galaxy, the best fitting R^/^ model, and the residuals. The residual 
images have been scaled by a factor of 10 compared to the data and models to bring out faint features. The 
bars at the bottom are one arcsecond in length. 



brightness profiles by fitting synthetic galaxy im- 
ages to our data using a range of surface bright- 
ness profiles, ellipticities and orientations. A se- 
ries of models were constructed using exponential 
surface brightness profiles, de Vaucouleurs R^/** 
profiles and the more general R^/" Sersic profiles. 
We did not consider more general fitting laws due 
to the relatively small range of radii (0'.'12 — 2", 
or 1 — 17 kpc at z ~ 1.5) covered by our ob- 
servations. Models with a range of scale lengths 
and eccentricities were convolved with the Point 
Spread Function (PSF) of the observations and 
subtracted from the NICMOS images. We used 
PSFs derived from well-detected unsaturated stars 
in each NIC3 field rather than the TinyTim sim- 
ulations as we found the former produced better 
fits. The residuals were computed and the model 
parameters were iterated to minimize the square 
of the residuals within the box of 8'.'4 x 8'.'4 cen- 
tered on each galaxy. The initial guess for the 
centroid was the position of the highest intensity 
pixel within the fitting box, and the total mag- 
nitude was estimated according to the total in- 
tensity confined in this box. Both initial guesses 
were made after masking out of the neighbouring 
sources. The root mean square (RMS) image was 
used to give relative weights to the background 
pixels during the fitting. By using different stars 
the width of the NIC3 PSF was allowed to vary 
to include the effects of spatial and temporal vari- 
ations in the NIC3 PSF. Changing the PSF had 
very little impact on the derived effective radii in 
all cases. The best-fit models for all galaxies in the 
sample are presented in Figure [1] (middle column) 
along with the residual images (last column). Pa- 
rameters of the best-fit R-'^/'* and R^/" profiles for 
each galaxy are given in Table [31 The listed min- 
ima of reduced x^ are well below unity, suggest- 
ing that the flux uncertainties introduced by the 
RMS images are overestimated. We performed 
the same morphological analysi s on the MUNICS 
galaxies (JLonghetti et al.ll2007t ). The NIC2 PSF 
used for modeling 2D profiles of these objects 
was derived from the TinyTim simulations. The 
resulting best-fit R^^^ profiles are graphically il- 
lustrated in Figure [21 The parameters obtained 
are listed in Table [H a long with the results from 
Longhetti et al.l (J2007I ) for comparison. The re- 
duced x^ are again below unity, but the values 
obtained for our best fit are very similar to the 



ones obtained for [Longhetti et al.l (|2007[ ) parame- 
ters, except for the total F160W magnitudes where 
the difference is greater then \a. The reasons for 
this discrepancy may be the simulated PSF we 
used for 2D fitting and the different methods ap- 
plied for background subtraction. Also, resulting 
R^/^ fit effective radius Re and surface brightness 
(/i)e for objects S2F1_142, S7F5_45, S2F1_633, and 
S2F1_443 differ for more then la from the previ- 
ously reported ones. When fitted with R-^/" pro- 
files, the best fits for the three of these objects - 
S2F1.142, S2F1.633, a nd S7F5-45 - have lowe r in- 
dices n than listed in [Longhetti et al.l ( 20071 ) - 2 
instead of 3.5, 2.5 instead of 4.1, and 1.5 instead 
of 2, respectively. On the other hand, the best 
fit Ri/" profile for S2F1_443 has index n = 2.8, 
higher than n = 1.9 reported bv [Longhetti et al. 
(|2007[ ). For the rest of the MUNICS sample the 



difference in th e goodness of fit f or R^/ ^ profile be- 
tween our and [Longhetti et al.[ (|2007l) analysis is 
A(x')<0.2. 

As a consistancy check, we also determined one- 
dimensional (ID) azimuthally averaged radial sur- 
face brightness profiles for each galaxy and for the 
corresponding models resulting from its 2D profile 
fits. These ID radial profiles were extracted using 
the approach developed by i Jedrzeiewskii (1987^ as 



implemented in IRAF (|Todv 1119931 ). Integrated 



magnitudes were determined within a series of el- 
liptical isophotes, the spacing of which grows with 
radius. We masked objects closer than 10 be- 
fore determining the surface brightness profiles of 
the galaxies. In most cases we are able to deter- 
mine the profile over roughly six magnitudes of 
surface brightness and to radii of 1'.'5, or ~ 13 kpc 
at z = 1.5. The 5cr limiting surface bright- 
ness for most of our observations is fipwo ~ 23 
mag arcsec"^; the data for 12-5869 and 12-5592 
reach approximately 0.3 magnitudes deeper. This 
surface brightness limit corresponds to fir ~ 20 
mag arcsec"^ (/i^ ~ 20.3 mag arcsec"^ for 12- 
5869 and 12-5592) for a galaxy at redshift z = 1.5 
that is formed at zy = 6 with exponentially declin- 
ing SFR and e- folding time r = 0.1 Gyr. Surface 
brightness profiles were determined in a similar 
fashion for each star that served as a local measure 
of the PSF. Azimuthally averaged surface bright- 
ness profiles for all of our GDDS objects are pre- 
sented in Figure [H with the profiles of best-fitting 
2D models and a PSF profile shown as solid lines 
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Fig. 3. — Upper panels: major axis surface brightness profiles in the F160W band for each galaxy (squares), 
with R^/* (green line) best-fit profile, R^'" (red line) best-fit profile, and a PSF profile (blue dashed line) 
overploted. The step used to present isophotal surface brightness corresponds to the pixel scale of our 
drizzled NfC3 images (0'.'12). The limiting surface brightness in each panel presents (roughly) 5cr limit for 
our observations. The lower part of each panel shows the residual differences between the data points and 
the model fits, with the 1 a errors on the data shown for comparison. 



and a dashed line, respectively. Figure [3] confirms 
that all galaxies in our GDDS sample are well re- 
solved, except for the target 12-6072 that seems 
only marginally resolved when compared to the 
PSF ID profile. The profiles are smooth in nearly 
all cases, the exception being object 15-4367 which 
shows a step at a = 1'.'5. Careful examination of 
this obejct's NIC3 image revealed that it was not 
perfectly symmetric and harboured a weak disk. 
The best R^/" profile index of ~ 2 confirms these 
findings. In addition, 15-4367 has a very faint 
neighbouring object that had to be masked out 
before fitting. These two effects produced the step 
in its ID profile seen in Fig. [31 

In order to estimate the errors on parameters 
obtained by our 2D and ID fitting procedures, we 
undertook a series of Monte Carlo (MC) simula- 
tions which incorporated all the sources of sys- 
tematic and random errors we were able to iden- 
tify. We constructed a set of galaxy images from 
our best-fit model for each galaxy and convolved 
these with a range of PSFs (i.e., PSFs obtained 
from different stars) and added these to the back- 
ground images. We dithered the position image 
about the central value to explore the importance 
of binning, and used RMS images to construct 2D 
arrays of random numbers to capture poisson noise 
and structure in the sky background. Each image 
constructed in this way went through the same 
fitting procedure as the real galaxy image from 
our sample. The standard deviations of resulting 
parameters are shown as the error estimates re- 
ported in Table [21 The reduced y^ values for the 
best fits to the MC simulations are of the order of 
unity and larger then reduced y^ of the best fits 
to the data, which makes our error estimates very 
conservative. 

3.2. K-corrections and cosmological dim- 
ming 

Our analysis requires comparison between the 
properties of our 1.2 < z < 2 samples observed at 
1.6 jiTu (i?-band) to those of present-day galaxies 
observed at visible wavelengths. In order to make 
a proper comparison, we need to transform the 
various data sets to a common bandpass and apply 
a K-correction. We computed appropriate spec- 
tral energy distributions (SEDs) using PEGASE- 
HR spectral synthesis models (JLe Borgne et al. 
20041 ). The model that we used is based on the 



JBaldrv fc Glazebrookl (2003) initial mass function 
(IMF), solar metallicity, and an exponentially de- 
clining star formation rate with a time scale of 
T = 0.1 Gyr, very similar to a single burst. The 
typical ages of GDDS and MUNICS passiv e galax- 
ies at 1.2 < z < 2 are 3 - 4 Gyr ( Paper IV . 
iLonghetti et al.ll2005[) and we used a 4 Gyr model 
to approximate their SED. It is important to em- 
phasize that the correction needed to reduce our 
iJ-band data to rest-frame SDSS-r is remarkably 
insensitive to SED shape since redshifted iJ-band 
closely matches rest-frame SDSS-r at z ^ 1.5. 
The photometry for the two samples is listed in 
Tables [31 and [H Cosmological surface brightness 
dimming will reduce the observed surface bright- 
ness and these must be corrected by {1 + z)"^ to 
transform them to the rest-frame. 

4. Results 

4.1. Morphologies of Passive Galaxies at 

z > 1.3 



All of the objects in our NICMOS F160W sam- 
ple (shown in Figures [H and [5]) have compact 
morphologies and none show obvious evidence of 
interactions, such as double nuclei or disturbed 
isophotes at bright levels. The star- forming mas- 
sive galaxies drawn from the GDDS sample, by 
contrast, exhibit a wi de range of disturbed mor- 
phologies as shown in [Paper Villi . The interme- 
diate age and composite population systems pri- 
marily have disk morphologies, w hile the pas sive 
galaxies at z < 1.3 discussed in [Paper Vllj ex- 
hibit a preponderance for compact and regular 
morphologies. Six of the 10 GDDS galaxies in the 
present sample appear to be early types with R^' " 
profile index n > 2.5 (Table [3]), while the four po- 
tential disk systems in our z > 1.3 passive sam- 
ple appear to have prominent bulges. Thus 60% 
of our GDDS sample defined by spectral proper- 
ties have pure early type morphologies, and this 
fraction rises to 90% when the prominent bulges 
with very weak disks are also taken into account 
as early type object. To a first approximation, our 
NICMOS Camera 3 images extend the correlation 
between spectra indicative of old stellar popula- 
tions and compact early-type morphologies from 
2;^1.3to2;~2. This is not surprising given 
previous indic ations in this direct ion from smaller 
samples fe.g.. [Cimatti et al.|[2004l) . 



Table 3 
Morphological parameters of the galaxies in the GDDS sample 
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Table 4 
Morphological parameters of the galaxies in the MUNICS sample 
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The correlation between color and morphologi- 
cal type is not as strong for the red galaxies, as 
a number of studies have shown. At redshifts 
near unity, red R — K or I — K selected samples 
contain roughly a s many disk as early-type galax- 
ies (e .g.jMousta kas et alj|2004 lYan fc Thompson 
20031 etc.). At higher redshifts red selected sam- 
ples also show a mix of morphologies, as shown 
fo r the z ^ 1.5 range inJPapcr VII j and at z > 2 



inJPai 

■ Istoc] 



bv lLabbe et"an(|2005[) . [stockton et an(|2004[ ). and 
others. 

4.2. Surface Brightness Profiles &; Sizes 

Azimuthally averaged surface brightness pro- 
files presented in Figure [3] confirm that six of our 
10 GDDS galaxies are well-fit by R^/"* profiles. 
The effective radii for these six objects range from 
as small as O'.'OS to as large as 0'.'42, or from 0.4 
to 3.6 kpc. The median effective radius is 0'.'26 
or 2.2 kpc. As Figure [T] shows, for the most part 
the 2-D models fit the data well and the resid- 
uals are not significantly greater than the sky 
noise. In 12-8895 and 12-5869 there appear to 
be some non-axisymmetric structures within the 
central one arcsecond, while in 12-6072 the model 
is too peaked. Four of our 10 GDDS galaxies are 
clearly better fit by R^/" profiles with indicies near 
2, rather than the R^/^ law. These are: 12-6072, 
12-8025, 15-4367 and 15-5005. As can be seen in 
Figure [3] the significance with which the R^/^ law 
fit is rejected in these objects is low except in the 
case of 12-8025 where the outer isophotes depart 
strongly from the R^/^ law profile. 

The effective radii of the GDDS galaxies are 
smaller than those of present-day cluster ellipti- 
cals and early- type field galaxies. The median ef- 
fective radius for low rcdshift cluster ellipticals is 
~ 4 kpc (jorgensen et al.lll995HSchombert Ill986[ l. 
and the fie ld early type galax ies at z ~ 0.5 from 
the CFRS ( Schade et al.ll999l ) have a fairly similar 
median size. The hosts of lu minous radio galax- 
ies at z ~ 0.8 - 1 studied by IZirm et all (|2003h 
probably represent the most massive end of the 
field & group early type populations at these red- 
shifts. Their sizes are also similar to the lower 
redshift samples and larger than the GDDS el- 
liptical galaxies that have median effective radius 
of 2.2 kpc. In contrast, the distant red galax- 
ies (DRGs), defined by their J — K colors, at 
2 < z < 3 have a median effective radius of 1.4 kpc 



( Toft et al.ll2007 ). somewhat smaller than the pas- 
sive GDDS galaxies in our sample at z ^ 1.7. 

The sizes of the GDDS passive galaxies ap- 
pear to support a fairly strong evolution in scale 
length among the early type galaxies in the 1 < 
z < 3 interval. A mundane potential explana- 
tion for this result is that the under-sampling of 
the NIC3 PSF data has led to imrehable fits. We 
can rule out this hypothesis on the basis of three 
tests. Firstly, we have re-fitted the six galax- 
ies with more finel y samp led NIG2 data from 
the iLonghetti et al.l ( 20071 ) sample, and we re- 
cover very similar fits (see Table |4]). These fits 
are shown in Figure |4] using das hed lines to join 
the va lues of points obtained by ILonghetti et al. 
(|2007l ) to those obtained by us. Secondly, we 



have undertaken detailed MC simulations (used to 
set our error bars in Figure [4|) based on generat- 
ing idealized over-sampled images which are ran- 
domly displaced by sub-pixel shifts before being 
binned to NIC3 resolution and re-fitted. Lastly, 
two of our objects - 12-5592 and 22-1983 - were 
observed in the F814W band with ACS on HST. 
The sizes that we measure for these galaxies, al- 
beit at shorter rest-frame wavelengths, are in good 
agreement with the sizes derived from our NIG3 
data. Thus we are confident that our size deter- 
minations are robust. 

The strong correation between mass and size, as 
measured by the effective radius, makes compar- 
isons between the average or median properties of 
different samples imprecise meas ures of evolution. 



The lower redshift samples (z < llJ(3rgensen et al 
19951 : ISchombert Ill986i : ISchade et al.lll999f l 



cover 



a broad range of the parent luminosity functions 
while the higher redshift objects (1 < z < 3), in- 
cludi ng the DRGs, the GDDS and MUNICS sam- 
ples (JToft et al.ll2007L IPaper Ilj . ILonghetti eFal 



20051 ) , sample the high mass end of the galaxy pop- 
ulation and thus are biased to large values in their 
median sizes. This further strengthens the conclu- 
sion that there is strong evolution in the charac- 
teristic sizes of early type galaxies above z ~ 1. 
The evolution in galaxy sizes can be further quan- 
tified by examining the size-mass correlation and 
its evolution, as is discussed in section [44l 



4.3. The Kormendy Relation to z = 2 

In Figure |4| we present the rest-frame r-band 
Kormendy relation, (/i)e vs. Re, for the GDDS 
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Fig. 4. — Mean rest frame Gunn-r surface brightness within effective radius Re as a function of Re (Kormendy 
relation) for objects at redshi fts 1.2 < z < 1.9 (GDDS and MUNICS samples) and for the sample of local 
galaxies (JBernardi et al.ll2003l SDSS,). The soHd Hne is the best-fit relation to the SDSS objects. The dashed 
lines represent the expected luminosity evolution of the local (SDSS, solid line) relation at z = 1.5 for galaxies 
formed at Zform = 6, 2 with exponentially declining SFR and e-folding time r = 0.1 Gyr. Different symbols 
correspond to different samples, and circled triangles denote re- fitted MUNICS sample. Left panel shows 
R^'"* profile parameters for the galaxies from the GDDS and MUNICS samples, while the right one shows 
their best-fit R^/" profile parameters. 



and MUNICS samples. As noted earlier, our 
construction of this diagram is particularly ro- 
bust because our observed 7J-band observations 
match rest- frame r-band at z — 1.5 and hence 
there is negligible residual K-correction uncer- 
tainty. We have not applied any evolutionary cor- 
rections to the observed surface brightness values. 
Figure |4] includes the corresponding distribution 
for present-day early- type galaxies from the SDSS 
(JBernardi et al.ll2003l ). 

Figure m shows that the tightness and slope of 
the Kormendy relation in the GDDS -I- MUNICS 
sample is similar to that defined by the local rela- 
tion. There is a hint that the high-redshift slope 
may be slightly steeper than the local value, but 
the difference is not significant. While the high- 
redshift ellipticals fall along a tight Kormendy re- 
lation, the relationship itself is offset to higher 
surface brightness from the low-redshift reference 
sample. The simplest explanation for this is that 
is that we are seeing galaxies nearer to their epoch 
of formation, when they are brighter, and thus the 
Kormendy relation is shifted upwards. This evo- 
lutionary effect cannot fully explain the evolution 
in the Kormendy relation. The offset in surface 
brightness compared to the z ^ sample is too 



large (^ 2.5 mag) to be explained by pure lumi- 
nosity evolution of stellar populations unless the 
redshift of formation is very recent (zform ^ 2), 
which is inconsistent with both their colors and 



spectra feee lPaper IVD which argue that these are 
old systems with Zform ^ 4. In the latter case, 
the maximum dimming allowed is 1 to 1.5 mag, 
depending on the selected IMF and the star for- 
mation history. In addition, we also see from this 
figure that, in spite of their large masses, typi- 
cal high-z ellipticals are substantially smaller than 
their local counterparts. In contrast to the median 
effective radius for the GDDS sample of 2.2 kpc, 
early-type galaxies in the SDSS sample presented 
in Fig. m span the range of effective radii with the 
median value of 4.9 kpc. Finally, we see that three 
out of ten galaxies in the GDDS sample are 'ultra- 
compact' (Re < 1 k pc), and thus are o f much 



higher stellar density. ICimatti et al.l pOOSh found 
a similar fraction from ACS imaging and estimate 
that the number density of comparably dense ob- 
jects at z = is up to 10'' times lower than at 
z — 1.5. In contrast, in the MUNICS sample of el- 
liptical galaxies (1.2 < z < 1.7) no 'ultra-compact' 
objects are found. As we will discuss in the follow- 
ing section, our findings lead us to also conclude 
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that strong size evolution (a factor of 2 or more) 
is the additional ingredient needed to explain the 
shift in the Kormendy relation. 

4.4. The Mass-Size Relation and the Ste- 
lar Mass Kormendy Relation 

As the previous section illustrates, a proper 
comparison between galaxy samples at high and 
low redshifts nearly always entails corrections for 
luminosity evolution. We can, however, improve 
on the standard procedure of using simple models 
of luminosity evolution by using multi-color SED 
data to fit stellar population models and derive 
stellar masses for the ga laxies in question (this 
was done and described in lPaper IIll for the GDDS 
sample). We then recast the data into a new 'stel- 
lar mass Kormendy relation' which allows a more 
fundamental comparison. By doing this we are 
using the complete set of information (the colors) 
to measure and remove the luminosity evolution. 
A further advantage to the use of stellar mass is 
that it allows us to compare optical and near-IR 
samples and plot them on the same diagram. A 
possible disadvantage is that we rely heavily on 
the mapping from light to stellar mass given by 
our spectral synthesis modeling, which, in turn, 
depends on the correctness of our assumptions. So 
for example derived masses would be in error if the 
assumed IMF is evolving rather than static. 

We consider two projections of the structural 
evolution that minimize the impact of luminosity 
and spectral evolution. The first is the size-mass 
relation, while the second is the relation between 
stellar mass density and size, which we will refer 
to as the stellar mass Kormendy relation. In de- 
riving the stellar mass density we assume that the 
F160W light traces the stellar mass. 

In Figure [5] we plot the size-mass relation 
for our sample. To enhance the usefulness of 
this figure, we augmented our GDDS and MU- 
NICS data using published measurements ob- 
tained for passive galaxies in the redshift range 
1.1 < z < 2.0 taken from from two surveys in 



the H UDF (jDaddi et al.l l2005bl : iMaraston et al. 



20061 ). a survey of six galaxies with dominant 
old stell ar population in the fie lds of radio-loud 
quasars (JMcGrath et al l |2007a||b|) , and GMASS 



(ICimatti et all l2008h . While (McGratheTZ 



20081) and HUDF (JDaddi et al.l l2005bl) effective 
radii were measured by fitting ACS F850LP (z 
band) galaxy images. We corrected all of the stel- 
lar mass determinatio n s to a common IMF, using 
Baldrv fc Glazebrookl (|2003|) IMF, a c cordi ng to 



the relati ons given in 



Cimatti et al.l (|2008l ) and 



Paper Ilj . Finally, to place our data in a broader 



context. Figure \5\ also shows the size-mass rela- 
tionship for local ear ly-type galaxies in the SDSS 



(jBernardi et al.ll2003D. We re c ompu ted the stellar 



masses for the lBernardi et al.l ( 20031 ) SDSS sample 
using th e same prescription applied to the GDDS 
sample (JBaldrv et al.l 120081 : IPaper IIll ). The de- 
ri ved masses are in good agreement with those 
of iKauffmann et al.l ( 2003! ) . The size-mass rela- 
tionship for early- type galaxies shown in Figured 
shows a number of interesting features, the most 
striking of which is that the high-redshift and low- 
redshift populations show relatively little overlap. 
In fact, they seem to describe nearly independent 
loci in size-mass parameter space, with similar 
slopes, but with galaxies at z ~ 1 — 2 systemat- 
ically smaller, at a fixed mass, than galaxies at 
z — 0. The error bars on individual data points 
are rather large, but taken as a whole, only ~ 25% 
of high redshift early-type galaxies lie in the re- 
gion of size-mass space occupied by low-redshift 
systems. 

The size-mass relationship of elliptical galax- 
ies at z ~ is well described by a power law 
with the same exponent (~ 0.5) as for the early- 
types at z ^ 1.5. Galaxies with stellar masses of 
8 X 10^° Mq, comparable to M* today, are approx- 
imately three times smaller at z ^ 1.5 than their 
apparent counterparts today. The number den- 
sity of compact galaxies with Rg < 1 kpc ('red 
nuggets') in the redshift range 1.1 < z < 2 is 
2 X 10~^ Mpc~'^. In contrast, number den sity of 



these objects in the SDSS sample (Bcrnard i et al. 
20031 ) is 3 X 10"^ Mpc^^, three orders of magni- 



([2007a) use NIC3 F160W observa tions for their 
morphological analysis, GMASS (jCimatti et al. 



tude lower than that for the higher redshift ob- 
jects. The 'red nuggets' in two samples are dif- 
ferent with respect to mass, too - the median of 
GDDS compact galaxies mass is 10'^"'^ M©, while 
the objects of the same compactness in the local 
Universe have masses with ten times lower median 
(i.e., 10^° M0). The passive galaxy population at 
1.1 < z < 2 span a similar range in stellar mass 
as galaxies today (2 x 10^° — 6 x 10^^ M©) so, at 
least at the high mass end, the bulk of the evo- 
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lution from z ^^ 2 to z r^ appears to be in size 
rather than mass. 

In Figure [6] we plot the projected stellar mass 
density within a radius equal to Rg (i.e, pe = 
3M*(R < Ro)/(47rR^)) versus Rg - the stellar mass 
Kormendy relation. This projection shows the 
evolution in the structural properties of the pas- 
sive early- type galaxies very clearly. The z > 1.1 
galaxies are offset to smaller radii and dramati- 
cally higher projected surface mass densities com- 
pared to massive early-type galaxies today. Com- 
pact objects in the local SDSS sample appear less 
dense since they are less massive than high red- 
shift objects with the same size. In the density 
space populated by red nuggets at higher redshifts 
{pe > 10^*^ Mq kpc"'^), there are no galaxies in 
the SDSS sample, implying that number density 
of these objects at z = is < 4 x 10~^ Mpc"'^. 

In both figures [5] and [6] we have color coded 
the symbols according to redshift into two sub- 
samples: 1.1 < z < 1.46 and 1.46 < z < 2. This 
splits the sample into two equal time intervals of 
duration 1.1 Gyr and nearly equal sample sizes. 
There is a signifcant diffence in the size distribu- 
tions in the two sub-samples. In the lower redshift 
sub-sample 6/18 galaxies, or ^ 33% of the sample, 
fall within the range of the local sample, while in 
the high redshift sample, only 4/25, or ^ 17% of 
the galaxies fall within the locus of the local sys- 
tems. Thus it appears that the strongest evolution 
in size is occuring in the 1 < z < 1.5 interval, al- 
though as we will describe in the next section, the 
heterogenous nature of the data does not allow us 
to conclude this with much confidence . A number 
of other studies (e.g.. lTreu et al.ll2005r i show that 
2 ~ 1 early type galaxies have normal sizes and 
mass densities. 

5. Discussion 

The key result of this paper are that the sizes 
and projected mass densities of early-type pas- 
sively evolving galaxies have changed very signif- 
icantly since z ~ 2. A number of other studies, 
noted above, have reached similar conclusions in 
samples with higher and overlaping redshift inter- 
vals. Our analysis has removed much of the un- 
certainty associated with evolutionary corrections 
in luminosity and spectral shape by dealing with 
the mass density rather than surface brightness. 



There are a number of potential explanations 
for the dramatic evolution in the sizes and densites 
of the passive galaxies. If the compact massive 
galaxies at z ~ 2 are to evolve into massive ellip- 
tical galaxies at z ~ they must grow by a factor 
of 2 — 3 in size. The two most plausible paths to 
this evolution are injection of energy into, or the 
loss of mass from, the central regions. One possi- 
bility is that mergers input energy into the stellar 
systems and increase their equilbrium sizes. The 
quiescient spectra of galaxies in the same stellar 
mass range at 1 < z < 1.5 suggest that any such 
merger be "dry" and produce little star formation 
and related activity. Dry mergers have been iden- 
tified as a likely evolutionary path for the com- 
pa ct massive galaxies a t z > 2 discussed recently 
by Ivan Dokkum et al.l ( 20081 ). The large stellar 
masses of the compact passive galaxies at z < 2 
suggest that equal mass mergers cannot be ubiq- 
uitous at later epochs. In Figures O and [S] we show 
vectors that approximate the impact of an equal 
m ass merger, based on the sim ulations performed 
bv iBovlan-Kolchin et al.l ( 20061 ). Galaxies become 
both larger and more massive and move primarily 
along the mass-radius and mass-Kormendy rela- 
tions rather than normal to them. This problem 
makes this explanation for size evolution unsat- 
isfactory. While there is good evidence for an 
increase of roughly a factor of two in the total 
stellar mass density in red sequence galaxies since 
z ~ 1.3, this appears to be in the form of new 
galaxies appearing on the red sequence rather than 
mass growth in previously passiv e systems (e.g., 
Faber et al.ll2007l : lBell et al.l2004[ ). One could per- 



haps appeal to many minor mergers to puff up a 
galaxy's size, but they would have to all be dry to 
keep a galaxy on the red sequence and numerous 
enough to have a significant effect, which seems 
somewhat contrived. 

It has been pointed out to us (N. Murray, pri- 
vate communication) that adiabatic expansion is 
an interesting alternative to dry merging for in- 
creasing the size of galaxies. This proces s has 
long been familiar to stellar dynamicists ( Hilla 
19801 ) and been verified by numerical simula- 



tion (e.g.. iBaumgardt et al.l 120071 ). The process 
has also been used to model the influence of 
strong stellar winds in conditio ning the G alac- 
tic globular cluster distribution (jZhad 120021 ). In 
the present context, the potential for adiabatic 
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Fig. 5. — Effective radius Ro as a function of stellar mass for five samples of early- type galaxies in the 
redshift range 1.1 < z < 2. Points are color-coded by two redshift ranges (red — z > 1.46, blue = z < 1.46). 
Different symbols correspond to different surveys, with triangles denoting re-fitted object from the MUNICS 
sample (as in Fig.Bt. The size-m ass relation for local early- type galaxies in the S PSS is presen t ed wi th sizes 
taken from iBernardi et al.l ( 20031 ). and matched with masses calculated following iBaldrv et al.l ( 20081 ) (black 
points). Contours represent linearly spaced regions of constant density of galaxies in size- mass parameter 
space. The solid line is the best- fit relation to the data poin ts at redshifts 1.1 < z < 2. Three arrows denote 
the effects that 1:1 dry merger ( Bovlan-Kolchin et al.l 120061 ). adiabatic expansion with 50% mass loss, and 
pure size evolution at constant stellar mass would have on the positions of both the least and the most 
massive galaxy. See text for details. 
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Fig. 6. — Stellar mass density within the effective radius Re as a function of Ro (the "stellar mass Kormendy 
relation") for five samples of early- type galaxies in the redshift range 1.1 < z < 2. Symbols are as in 
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point s at redshifts 1.1 < z < 2. Three arrows denote the effects that 1:1 dry merger (jBovlan-Kolchin et al. 
20061 ). adiabatic expansion with 50% mass loss, and pure size evolution at constant stellar mass. See text 
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expansion to explain the existence of massive 
small ellipticals at higli redshift is developed 
in a paper by iMurrav. Quataert fc Thompsonl 



(|2008. in preparation! , hereafter MQT). To mo- 
tivate the present discussion, a basic version of 
the some of the key theoretical ideas in the latter 
paper, kindly communicated to us in advance of 
publication by the authors, will be applied to the 
GDDS sample here. 

Adiabatic expansion will occur in any relaxed 
system that is losing mass. As mass is lost the po- 
tential becomes shallower, so the system expands 
in order to relax into a new stable equilibrium. 
The amount that a system expands depends on 
both the e xtent and speed of the mass loss (see 
Zhadl2002l for details). In general, if a fraction 
^ = (^initial - "ifinaO/n^initiai of the total mass 
is lost on a dynamical timescale (or longer), the 
size of the system increases by a factor of approx- 
imately Am ■ If the mass is lost more quickly 
than the dynamical timescale, then the expansion 
of the system will be larger than this estimate. It 
is trivial to show that as the system loses mass 
the dynamical timescale increases in proportion 
while the escape velocity decreases as 



to 



(1-^) 



1 — ^21 so there are at least two sources of posi- 
five feedback leading to further increase the size as 
the system evolves. Of course, in the extreme case 
where a significant fraction of the total mass is 
lost on a short timescale, the system may become 
unbound. 

What processes might lead to mass loss in el- 
liptical galaxies? The obvious candidate is stellar 
winds from sites of active star formation. How- 
ever, the early-type galaxies being studied here 
are relatively red and spectroscopically passive, 
so winds from young stellar populations are un- 
likely candidates for mass loss. An interesting al- 
ternative is mass loss from evolved A and F-type 
stars, and we have explored this ideas using the 
following toy model. We model a galaxy as an 
instantaneous burst with a solar-metallicity stel- 
lar population whose main sequence lifetime (as 
a function of mass) is that given in Table 5.2 of 
Binney & Merrifield (1998). We assume that after 
leaving the main sequence all stars more massive 
than 8 solar masses wind up as stellar remnants of 
1.5 solar mass, and that all stars less than 8 solar 
masses wind up as remnants with 0.6 solar mass. 



We also assume that mass loss from stars is never 
recycled into future star formation and it onflows 
far out into the galaxy's potential well, or is lost 
completely. 

In this case, -^^ as a function of time takes 
on the form shown in Figure [7] for three initial 
mass functions (Salpeter IMF, Scalo IMF, and the 
Baldry & Glazebrook (2003) IMF). Our toy model 
suggests that -^^i rises sharply with time until 



ages of around 2 Gyr, at which point 



Am(t) 



flat- 



tens out, peaking at around 30% for the Salpeter 
IMF, and at 50% for the top-heavy Baldry & 
Glazebrook (2003) IMF. Thus the degree of mass 
loss from a very top-heavy IMF could explain the 
size growth. This is shown by the arrows in Fig- 
ures [5] andj6l_which_sliow_tlie efects that 1:1 dry 
merger ( Bovlan-Kolchin et al.ll2006l . cyan arrow), 
adiabatic expansion with 50% mass loss (magenta 
arrow), and pure size evolution at constant stel- 
lar mass (green arrow) have on the positions of 
both the least and the most massive galaxies in 
our sample. However, the timescale over which 
this occurs poses a huge challenge for explain- 
ing the size growth entirely by adiabatic expan- 
sion. In this paper we study the size distribu- 
tion of the population at a time when their the 
stellar populations are already rather old (see Ta- 
bles [T] and [2] and discussion in IPaper IVl ) so over 
the redshift range being probed the galaxies are 
old enough that the mass loss curves in Figure [7] 
are already nearly flat. Another constraint on the 
importance of adiabatic expansion is that is does 
not explain the steady factor of (at least) three 
growth in the stellar mass density locked up in 
massive gal axies over the redshift range 1 < z < 2 



reported in iPaper IVl and in other surveys , (e.g., 
iDickinson et all 120031 : iRudn ck et aD 120061) . espe- 
cially on the red sequence ( Paper VII]] ). As the 
typical mass does not appear to evolve (Fig. [5|) 
this primarily seems to be an evolution in num- 
ber. 

In spite of the problems noted above, adiabatic 
expansion does appear attractive because it moves 
the high-redshift distribution shown in Figure [5] in 
the right direction to match the low-rcdshift dis- 
tribution shown in the figure. This is not the case 
with equal-mass dry mergers, which, as shown by 
the cyan arrow s in the figure, and a s noted by pre- 
vious authors (JBovlan-Kolchin et a l. 2006), drive 
evolution along the Kormendy relation rather than 
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Fig. 7. — The mass loss fraction ^^^i as a function of population age in Gyr, for the simple model described 
in the text. We assume an instantaneous burst of star formation and show ^^ as a function of time with 
three initial mass functions. As expected, the total mass lost is a strong function of the fraction of stars at 
the high mass end of the IMF. The relative mass loss is small in the age range 1 — 7 Gyr (dashed lines). See 
text for details. 



displacing the relation itself. While a top-heavy 
IMF loses enough mass to grow the galaxies by 
the required factor of two over their complete life- 
time, the main problem with the adiabatic ex- 
pansion model is that to explain our observations 
that mass loss would have to occur over the age 
range of 1 — 7 Gyr, over which Fig. [7| shows only 
a 5 — 10 % effect. Ages of the GDDS galaxies 
are taken from Paper IV, and it is worthwhile to 
consider whether we might have significantly over- 
estimated the ages of the galaxies in that paper. 
Wc think this unlikely for two reasons. Firstly, be- 
cause broad-band color-based ages for these galax- 
ies seem consistent with ages inferred from spectra 
of these systems, which often exhibit photospheric 
features from old stars. Secondly, because chang- 
ing to a mor e top-heav y IMF than the Saipeter 
IMF used in IPaper IVl would not result in sys- 
tematically younger ages. In fact the reverse is 
true, since a more top-heavy IMF would tend to 
produce synthetic spectra which are bluer for a 
given star formation history at a given age. So 
to match the observed colors, any fitting routine 
would compensate by deriving older ages for the 
best fit. Quantitatively, we checked the size of 
this effect by generating models with an exponen- 
tially declining star-formation history (e-folding 
timescale t = 1 Gyr) with various stellar metallici- 
ties, using both Saipeter and BG03 IMFs (without 
extinction). We determined that ages using the 
(top-heavy) BG03 IMF are ~ 40 - 50% larger for 



galaxies which are found to be ~ 1 Gyr old using 
a Saipeter IMF. (Note that derived metallicities 
using the BG03 IMF are larger too). 

Some constraints on the duty cycle for the size 
change can be inferred from our observations, by 
noting that the redshift range spanned by our sam- 
ple is 1.1 < z < 2.0, corresponding to a spread in 
time of ^ 2.2 Gyr. The division of the sample in 
half at z — 1.46 using different symbols in Fig- 
ures [5] and [6] subdivides this redshift interval into 
two equal time bins, each of which is ^ 1.1 Gyr 
wide. The sample shown in Figures [5] and [6] con- 
tains data from a number of different surveys, and 
it is certainly unwise to attempt to compare the 
high-redshift and low-redshift subsets at a detailed 
level. But it is perhaps worth noting the following 
very general qualitative trends. Figure [5] appears 
to show that the character of the size-mass dis- 
tribution is rather different in the 1.1 < z < 1.46 
and 1.46 < z < 2.0 intervals, with neither distribu- 
tion resembling the local data distribution closely. 
This suggests some degree of evolution between 
the bins, but with the caveat that these two red- 
shift bins primarily consist of data from different 
surveys so the strength of the evolution cannot be 
confidently inferred. On a more speculative note, 
it can be argued that nothing in Figure [5] rules 
out the possibility that the high-redshift distribu- 
tion is evolving into the low-redshift distribution 
differentially, with different physics operating at 
the low mass and high mass ends. In fact, some 
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evidence for this is also hinted at in the figure, 
which appears to show that the smallest and least 
massive galaxies lie at 2; > 1.5. It is possible that 
dry mergers may well be growing the smallest and 
least massive galaxies along the fundamental plane 
early in a galaxy's life cycle, before some other pro- 
cess takes over and grows them further in some 
other way. 

It is interesting to contrast the data pre- 
sented in Figures [5] and [H] with data which spans 
the redshift range in between the SDSS data 
and our high-redshift observations. Figures [8] 
and [9] augment the data in Figures [5] and [6] 
with intermediate-redshift data taken from the 
CFRS dSchade et al. 19991: iLillv et all [l995) and 
GOO DS /DEIMOS (fBimdv et al.ll2007HTreu et al. 
120051 ) surveys. Effective radii for the CFRS ob- 
jects are obtained from the WFPC2 814W im- 
ages. Estimates based on the images in three 
ACS filters (606W, 814W, and 850W) are avail- 
able for the GOODS/DEIMOS sample. All ob- 
jects shown in the upper three panels of figures [8] 
and M have sizes based on the WFPC2 or ACS 
814W imaging that translates approximately into 
the rest-frame T^-band for the median redshifts in 
the 0.2 < z < 0.5 and 0.5 < z < 0.7 panels, and 
into the rest-frame i?-band for the median redshift 
in the 0.7 < z < 1. panel. GOODS/DEIMOS ob- 
jects in the 1 < z < 1.46 panel are presented with 
the effective radii in ACS 850W filter (approx- 
imately i?-band rest frame). The CFRS masses 
are obtained foUowing lBaldrv et al.l (2008j and us- 
ing imaging data of relatively low quality. We note 
that the difference in the rest-frame wavelengths 
that are probed at different redshifts makes it 
impssible to draw any quantitative conclusions 
about galaxy size evolution. However, figures [8] 
and [9] show qualitative trends consistent with 
smooth evolution over the 0.2 < z < 2 range. The 
dispersion on the size-mass plot in the 0.2 < z < 1 
regime is large (upper panels in Figures [8] and [9]) , 
but there seems to be some evidence for a sys- 
tematic offest relative to the local trends with the 
increasing redshift. The GOODS/DEIMOS data 
in Figures [8] and [9] spans both the low-redshift and 
high-redshift loci identified in each panel of Fig- 
ures |5] and [H] by contours and the line of the best 
fit, respectively. However, the majority of the 
low-redshift (0.2 < z < 0.5) GOODS/DEIMOS 
data lie closer to the local relation, in contrast 



to the 0.7 < z < 1 panel where the most of the 
GOODS/DEIMOS points are close to the z - 1.5 
objects locus. The CFRS data in the 0.7 < z < 1 
panel of Figure [5] does not seem to follow this 
trend, and we suggest that it may be due to 



the s hallow imaging of these objects (jLillv et al. 



19951 ). In the lower redshift panels of Figure [8] 



(z < 0.7) the positions of the CFRS objects are 
consistent with the GOODS/DEIMOS dataset. 
In order to compare the number of high mass ob- 
jects at different redshifts, we use a subsample of 
68 GOODS/DEIMOS objects with masses above 
the GDDS detection limit (see section 12. ip . It 
is interesting to note that relatively few (14/68, 
~ 21%) points from the GOODS/DEIMOS sub- 
sample have masse s greater than 1.5 x 10^^ Mq 
(M* at 1 < z < 2. lFontanaet al.ll2006l ). In con- 
trast, the high-redshift data set presented in Fig- 
ures [5] and [6] includes large fraction of objects 
with M^ > M* - 18/43 (- 42%). While this 
could perhaps be consistent with adiabatic mass 
loss, the arguments presented in our discussion 
of Figure [7] are compounded by the data pre- 
sented in Figures [8] and [9] which indicates that 
size growth is still occurring in galaxies even older 
than those in our GDDS sample. We think it is 
likely that the absence of very high mass objects 
in the GOODS/DEIMOS data is simply due sub- 
tle differences in various groups' methodologies for 
computing stellar masses from photometric data. 
To further address the question of structural evo- 
lution of galaxies presented in Figures [8] and [9] 
we plot the redshift dependence of the projected 
stellar mass density (defined in section [44|) in Fig- 
ure [lOl Dashed lines encompass the range of mass 
density which contains 90% of the local (SDSS) 
data points. The median stellar mass density of 
the SDSS galaxies is pe = 1.1 x lO^Mg kpc'^ 
and this value is denoted by a red cross plotted 
at z = 0.1 in Figure [TOl Large fraction (88%) 
of the GOODS/DEIMOS objects have mass den- 
sities above the local median value, and 65% of 
these galaxies have mass densities above the up- 
per dashed line in the figure. For the 1.1 < z < 2 
sample the corresponding numbers are 90% and 
77%, respectively. On this basis we can conclude 
that the stellar mass density increases over an ex- 
tended redshift range, though the dispersion of the 
plot is large, and more points in both intermediate 
and high redshift regime are needed to properly 
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constrain this redshift dependence. We intend to 
revisit the topic in a future paper. 

On balance, we conclude that at present neither 
adiabatic expansion nor equal-mass dry mergers 
seem able to explain the size growth in early-type 
galaxies. A successful model will have to simul- 
taneously explain the size change in the galaxies, 
the duty cycle for this size change, and the epoch 
in a galaxy's life history at which the change oc- 
curs. And, as noted above, mass density growth 
over the redshift interval being probed suggests 
that the size growth being witnessed is operat- 
ing within a broader context for galaxy formation. 
Over the redshift interval where early-type galax- 
ies are growing in size, the volume-averaged stellar 
mass density in massive galaxies is increasing, and 
the morphological mix is changing. 

6. Conclusions 

The size-mass relationship for early-type galax- 
ies evolves significantly from z = 2 to z = 1. Over 
the whole of this redshift range early type galax- 
ies tend to be a factor of 2 — 3 smaller than local 
counterparts of similar mass. Similarly compact 
galax ies are seen at z > 2 (jvan Dokkum et al. 



20081 ) ■ and we speculate that the very compact 
galaxies studied in the present paper are simply 
the evolved counterparts of these higher-redshift 
objects, caught at a time before subsequent size 
growth. By comparing the size distribution of 
our sample with that of lower redshift surveys, we 
conclude that significant size growth is probably 
occurring over the redshift range explored in the 
present paper. The physics of this growth remains 
mysterious. By comparing the size-mass relation 
at z '^ 1.5 with its local counterpart we conclude 
that equal mass dry mergers play only a limited 
role in growing early-type galaxies, at least once 
they are older than a few Gyr. Other processes 
may be as important as dry merging in growing 
early-type galaxies. Adiabatic expansion is one 
such process that we have examined, and while 
it may be important in growing young early-type 
galaxies, it is hard to sec how this mechanism can 
be invoked to obtain a factor of two growth in the 
sizes of galaxies as old as those in the present sur- 
vey. 
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